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' The axion is the most attractive candidate to solve the strong CP problem in 
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' QCD. If it exists, the inflationary universe produces axion fluctuations which 
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are mixtures of isocurvature and adiabatic fluctuations in general. We investi- 
gate how large isocurvature fluctuations are allowed or favored in order to ex- 
plain observations of the large scale structure of the present universe. Generic 
flat universe models with mixed (isocurvature-|- adiabatic) density fluctuations 
are studied. It is found that the observations are consistent with the mixed 
fluctuation model if the ratio a of the power spectrum of isocurvature fluc- 



^ ' tuations to that of adiabatic fluctuations is less than ~ 0.1. In particular, 



the mixed fluctuation model with a ~ 0.05, total matter density VLq = 0.4, 
and Hubble parameter Hq = 70km/s/Mpc gives a very good fit to the ob- 
servational data. Since the height of the acoustic peak in the angular power 
spectrum of the cosmic microwave background (CMB) radiation decreases 
significantly when the isocurvature fluctuations are present, the mixed fluc- 
tuation model can be tested in future satellite experiments. Ratios of the 
amplitude at the peak location to that at the COBE normalization scale 
for various models are given. Furthermore, we also obtain the amplitude of 
isocurvature fluctuations as a function of axion parameters and the Hubble 



parameter during the inflation. We discuss the axion fluctuations in some re- 
alistic inflation models and find that a significant amount of the isocurvature 
fluctuations are naturally produced. 
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I. INTRODUCTION 



The axion is a Nambu-Goldstone boson associated with spontaneous breaking of the 
Peccei-Quinn U{1) symmetry which is introduced in order to solve the strong CP problem 
in quantum chromodynamics (QCD). The breaking scale Fa of the Peccei-Quinn symmetry 
is stringently constrained by accelerator experiments, astrophysics and cosmology, and it 
should lie in the rage 10^'^ — lO^^GeV in standard cosmology. In particular, the upper limit 
is set by requiring that the cosmic density of the axion does not exceed the critical density 
of the universe. In other words, the axion should represent the dark matter of the universe 
if Fa is ~ lO^^GeV, which makes the axion very attractive. 

In the inflationary universe, quantum fluctuations of the infiaton field result in adiabatic 
density fluctuations with a scale-invariant power spectrum which would account for the large 
scale structure of the present universe. This natural generation of density fluctuations is 
one of successes of the inflationary universe. However, if the axion exists, another kind of 
fluctuations is produced in the inflationary universe. During inflation the axion has quantum 
fluctuations whose root mean square amplitude 6a is given by H/(2'k). {H is the Hubble 
parameter during inflation.) The quantum fluctuations of the axion are stretched by the 
inflation and become classical. When the axion acquires a mass ma at the QCD scale, these 
axion fluctuations result in density fluctuations of the axion. However, since the axion is 
massless during inflation, axion fluctuations with wavelength larger than the horizon size do 
not contribute to the total density fluctuations. For this reason, such density fluctuations are 
called "isocurvature" . Therefore, if indeed the dark matter consists of axions, it has both 
adiabatic and isocurvature fluctuations and may play an important role in the structure 
formation of the universe. 

In a previous work the large scale structure formation due to the mixed (adiabatic -|- 
isocurvature) fluctuations of the axion was studied with matter density parameter fio = 1. 
It was found that the introduction of isocurvature fluctuations significantly reduces the 
amplitude of the power spectrum P{k) of the density fluctuations, keeping P{k) normalized 
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by the COBE data. Since it is known that the cold dark matter (CDM) model with pure 
adiabatic density fluctuations and fio = 1 predicts an amplitude of P{k) that is too large on 
galaxy and cluster of galaxies scales, the mixed fluctuation model gives a better fit to the 
observations, although the shape P{k) on larger scales does not quite fit the observations. 
Furthermore, it was also pointed out that the height of the acoustic peak in the angular 
power spectrum of the cosmic microwave background (CMB) radiation decreases when the 
isocurvature fluctuations are added, which can be tested in the future satellite experiments. 
Similar studies were also done in Refs. 2) and 3).[] 

However, previous works deal with only restricted cosmological models (i.e. fio = 1). It 
is well known that the shape of the power spectrum inferred from the galaxy survey favors 
a low matter density universe. Therefore, in this paper we consider the structure formation 
with the mixed fluctuations in a generic flat universe (i.e. f2o + Aq = 1, where Aq is the 
density parameter for the cosmological constant) and investigate how large isocurvature 
fluctuations are allowed or favored in order to explain the observations. The constraint on 
the amplitude of the isocurvature fluctuations is reinterpreted as a constraint on the Hubble 
parameter H and the Peccei-Quinn scale Fa- We also study the cosmological evolution of 
the axion fluctuations in some inflation models. Among many inflation models, we adopt a 
hybrid inflation model and a new inflation model. In the hybrid inflation model, which is 
most natural in supergravity, we find that a significant amount of isocurvature fluctuations 
of the axion are naturally produced if Fa — 10^^~^^ GeV. We note that late-time entropy 
production may easily raise the Fa up to ~ lO^^GeV. In the new inflation model, we also find 
that isocurvature fluctuations are produced with a canonically favored value of Fa ~ lO^^GeV 
without late-time entropy production. 

In §|T| we investigate the structure formation with mixed fluctuations and compare the 
theoretical predictions with the present observations. The generation of isocurvature and 



A similar result was obtained in the context of the hot dark matter model. |^ 
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adiabatic fluctuations of tlie axion is described in §|TT|. We discuss some inflation models and 
isocurvature fluctuations of the axion in § [rV| . §|V| is devoted to conclusions and discussion. 
Throughout this paper, we set the gravitational scale ~ 2.4 x lO^^GeV equal to unity. 

II. COMPARISON WITH OBSERVATIONS 

The density field of the universe is often described in terms of the density contrast, 
6{x) = 6p{x)/p = {p{x) — p) /p, and its Fourier components Sk = y J (fx S{x) exp (ik ■ x) , 
where p is the average density of the universe, x denotes comoving coordinates, and is a 
sufficiently large volume. If 5{x) is a random Gaussian field as predicted by the inflation, 
then the statistical properties of the cosmic density field are completely contained in the 
matter power spectrum P{k) = V{\6k\'^), where (■ ■ ■) represents ensemble average. 

In general, the inflation predicts an adiabatic primordial spectrum, -Pad*™(^) ^ ^"'"■> 
where k is the wavenumber of perturbation modes, and a spectral index. In this paper, 
we set Hg = 1 (the Harrison-Zeldovich spectrum), which is predicted in a large class of 
inflation models including hybrid inflation, new inflation, and the chaotic inflation models 



discussed in § [rV| . As for isocurvature perturbations, on the other hand, the primordial 



spectrum index is conventionally expressed in terms of entropy perturbations Sar = 5a — |5r 
as Sar{k) oc where 6a and 6r are density fluctuations of axion and radiation fields, 
respectively. Employing this definition, we can express the primordial power spectrum as 
Piso^{k) oc k"'"'^^. We set = —3, which corresponds to the Harrison-Zeldovich spectrum. 

The primordial power spectrum is modified through its evolution in the expanding uni- 
verse and the present power spectrum can be described as 

Pik) = Pad + Piso = A^^kT^^ik) + A,,,kTl,{k), (1) 

where A^^ and Ai^o are normalizations of adiabatic and isocurvature perturbations, and 
Ta,d{k) and T[so{k) are transfer functions of adiabatic and isocurvature perturbations. Con- 
ventionally, the transfer function of the isocurvature perturbations is defined for the initial 
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entropy perturbations. Therefore the present matter power spectrum may be written as 
-Piso = Asok'^T^oif^)- Fo^^ this definition, hmfc_oo 7iso(^) = 1, while \imk^oTiso{k) = 1 for 
Eq.(|I]). However, this definition makes direct comparison between adiabatic and isocurva- 
ture matter power spectra difficult. Therefore we employ Eq.(|I]) to define the amplitude and 
the transfer function for both adiabatic and isocurvature perturbations. 

Here we define the ratio a of the isocurvature to the adiabatic power spectrum as ||l|J 



A. 

^iso /„\ 

a = — - . (2) 

It is well known that the transfer functions for adiabatic CDM models are essen- 
tially controlled by a single parameter, F = QqH, or more precisely, F = 
^0^(^0/2.7^)"^ exp(— fi^ — \^2hilB/^o)- m Here h is the present Hubble parameter nor- 
malized by 100km s~^ Mpc~^, To is the present cosmic temperature, and fi^ is the ratio 
of the present baryon density to the critical density. From the galaxy survey Peacock and 
Dodds [| estimated that F ~ 0.25 ± 0.05 + 0.32(n;i - 1). 

In observational cosmology, the quantity as, which is the linearly extrapolated rms of 
the density field in spheres of radius 8/i"^Mpc, is often used to evaluate amplitudes of the 
density fiuctuations. This is motivated by the fact that the rms fiuctuation in the number 
density of bright galaxies measured in a sphere of radius of 8/i^^Mpc is almost unity. 
Employing the top hat window function, we obtain 

dk ^ ( 3ji{kro)\ 



(3) 

ro=8/i-lMpc 



2-7?^ J k \ kvo 

where ji is the spherical Bessel function. 

Anisotropics of the cosmic microwave background radiation (CMB), 6T/T, can be ex- 
panded as 

^(7) = E E (^lmYim{l), (4) 
£=2m=-£ 



^This parameter can be expressed by using the Hubble parameter during inflation and axion 
breaking scale, as Eq.(|^). 
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where Y^m is a spherical harmonic function, and 7 denotes the direction in the sky. The tem- 
perature autocorrelation function (which compares the temperatures at two different points 
in the sky separated by an angle 6) is defined as {^^(7)^(7')^ = ^)CeP£{cos6), 
where 7- 7' = cos^, and Pi is the Legendre polynomial of degree i. The coefficients Ce are 
the multipole moments: 

e 

{2i+l)Ce={ E M'). (5) 
m=-e 

The predictions for the CMB anisotropics can be obtained by numerical integration of the 
general relativistic Boltzmann equations. 

Since the isocurvature fluctuations are independent of adiabatic fluctuations and give 
different contributions to the CMB anisotropics, the CMB multipoles can be expressed as a 
linear combination of two components as 

C,(«) = g{a)Cf + h{a)Cr , (6) 

where C^'^ and Cf° are adiabatic and isocurvature components of CMB multipoles, respec- 
tively, and g{a) and h{a) are functions of a. Note that C^'^ = (7^(0) and = Cf(oo). 
Employing the COBE normalization obtained from the 4-yr data, 0] which fixes Cio, i.e., 
Cio = C'lo = C'{q°, we obtain h{a) = 1 — g{a). Then we introduce a relative ratio between 
the isocurvature and adiabatic components at a = 1 as ff = (1 — g{l)) Cf"/g{l)Cf^. In the 
large scale limit (£ — > 2), this factor is almost fi ^ 6, which is expected from the differences 
between Sachs- Wolfe contributions of adiabatic and isocurvature fluctuations. |Tl| On the 
COBE scale ~ 10), however, is smaller than 6. Furthermore, for larger £, strongly 
depends on I because the angular power spectra for the adiabatic and the isocurvature 
fluctuations differ greatly (see Fig.||). The precise value of /lo depends on VLq and h. For 
example, /fg — 30 for fio = 1 and h ^ 0.6, and it becomes smaller when VLq or h becomes 
smaller. Using this factor Z^, we can write the ratio of the isocurvature to the adiabatic 
mode at arbitrary a as 

, (1 -g(a)) cr 
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Using the COBE normalization, we get g{a) = 1/(1 + fiQ<y). Eventually, if we know the 
value of /lo, the COBE normalized pure adiabatic (C^*^) perturbations, and isocurvature 
(Cf °) perturbations, we can obtain the COBE normalization for any admixture of adiabatic 
and isocurvature perturbations as 



= +/if . (8) 



A similar argument can be applied to (Xg as 



(Ts[a) 



\ 1 + /fo« 

where a^'^ and crg° are the values of ag normalized to COBE for the cases of pure adiabatic 
and pure isocurvature perturbations, respectively. 

We can obtain erg as a function of Qq from observations of the cluster abundance, since 
this abundance is very sensitive to the amplitude of the density fluctuations. Here we adopt 
the values of ag which are obtained from the analysis of the local cluster X-ray temperature 
function: |12| 



as = (0.52 ± 0.04)^0 (for ^]o + Aq = 1). (10) 

In other analyses, [[l^ -[l5|| similar values for erg have been obtained. 

In Figs.|] and |^, we plot erg (a) as functions of Qq and h. The observed Cg mentioned 
above is also drawn as a region inside the two short dash - long dash lines. Furthermore, 
three dashed lines, corresponding to the shape parameters F = 0.2, 0.25, and 0.3 from left to 
right, are drawn in the flgures. From these flgures, one can see that, for example, if /i ~ 0.7 
and F = 0.25 ± 0.05, the observed ag is inconsistent with the COBE normalized predictions 
for a = 0. It is possible, of course, to flt the observation with a = if /i ~ OA^Qq ~ 0.6). 
However, small values of the Hubble parameter {h ^ 0.4) seem unlikely from the recent 



observations, [0 which suggest h = 0.73 ± 0.06(stat) ± 0.08(sys). Rather, a = 0.04 to 
0.06 give better flts for h ~ 0.7. This fact is also seen in Fig.^ (left panel), where the 
power spectrum for a = 0.05, Qq = 0.4 and h = 0.7 is shown together with observational 
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data. In Fig.^ we also plot the pure adiabatic power spectrum for Qq = 0.3 and h = 0.7 
for comparison (right panel). The values of Qq for both power spectra are chosen to satisfy 
Eq. (|T0|) for h = 0.7. It is clear that the mixed model {a ^ 0.05) gives a much better fit to 
the data obtained from the galaxy survey. 

For each Qq and h, we can obtain the best fit values of a, which are plotted in Fig.^. The 
three dashed lines correspond to F = 0.2, 0.25, and 0.3, from left to right. One can see from 
this figure, for example, for h = 0.7 and F = 0.25, the best fit value of a to observations is 
a ~ 0.05. This figure also indicates that if F = 0.25 ± 0.05, the observations of ag require 
a ^ 0.1. This is also seen in Fig.||. 

In the above discussion, we see that a ~ (9(10~^) is consistent with observations of large 
scale structures. Let us examine what this value predicts. As mentioned earlier, the isocur- 
vature fluctuations give different contributions to the CMB anisotropics than the adiabatic 
fluctuations. If one takes very large values of a, the prominent acoustic peaks disappear, 
and only the Sachs- Wolfe plateau exists in the angular power spectrum. Therefore, we must 
investigate the height of the acoustic peak(s) predicted by our model. 

In order to predict this height quantitatively, we introduce the peak-height parameter in 
the CMB anisotropy angular power spectrum defined by 

i{i + l)Ce\,^^o ' ^ ^ 

where 7, Ci and ipeak are all functions of a. Here i{i + '^)Ce\i^i^^^^ denotes the highest value, 
which is usually the value at the first acoustic peak. When a becomes larger, the height of 
the acoustic peak becomes lower. If the peak height becomes lower than the Sachs- Wolfe 
plateau, this peak-height parameter tends to be unity. In Fig.H, we plot this quantity. One 
can see that, for example, for a ~ 0.05, = 0.4, and h = 0.7, this ratio is about 7 ~ 2 - 3. 

In Fig.|^ we plot the CMB angular power spectra normalized by the COBE for the cases 
a = 0, 0.05, and 00. Also, in Fig.||, we plot the same angular power spectra, changing the 
cosmological parameters Qq and h. From these figures, one can see how much the height of 
the acoustic peak decreases when a is large. 
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Future satellite experiments such as MAP and PLANCK [|T8| are expected to measure 
the CMB anisotropies with fine angular resolutions and detect the peak-height parameter. 
If this is done, our model can be tested by these observations. 



III. FLUCTUATIONS OF AXIONS 

During the infiation, the axion experiences quantum fiuctuations whose amplitude is 
given by 

iHk)f = ^, (12) 

where k is the comoving wavenumber and H the Hubble parameter during the infiation. 
These fiuctuations become classical due to the exponential expansion of the universe. After 
the axion acquires a mass at the QCD scale, the axion fiuctuations lead to density fiuc- 
tuations. The density fiuctuations of axions are of an isocurvature type, because they are 
massless during the infiation and hence do not contribute to the fiuctuations of the total 
density of the universe. The density of axions is written as 

Pa = \mla^ = \mlFlel (13) 

where is the mass of the axion and 9a the phase of the Peccei-Quinn scalar (0 < 0^ < 2ti). 
Then the isocurvature fiuctuations with comoving wavenumber k are given by 

cw-f^w) = (14) 

\Pa J a ta^a 

where we redefine H as the Hubble parameter when the comoving wavelength k~^ becomes 
equal to the Hubble radius during the infiation epoch. It should be noted that the above 
initial spectrum is of the Harrison- Zeldovich type in the case of isocurvature perturbations 
(see §11). 

On the other hand, the infiaton itself generates adiabatic fiuctuations given by 
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where V is the potential for an inflaton, and R{t) and H are the scale factor and Hubble 



constant at some arbitrary time t. [|T9| To compare these two types of fluctuations, it may 
be natural to consider the ratio of the power spectra at horizon crossing, i.e. k~^R = H~^, 
which is written as 



Pi. 



ISO 

ttKSY 



f\2 



9{V , 

„n=H ^ ^H^F^er ^^^^ 



Pad 

However, this definition is different from the ratio a of the present power spectra in Eq. 
since we need to take into account the difference between the time evolution of adiabatic 
and isocurvature perturbations for the radiation dominant regime and the matter dominant 
epoch. [|TT| Eventually, we can obtain the following relation: 00 

2\2/io\2 /4x2 4(y/^2 

rtTt-QV = CkT^QV = — 



" = [-5) IyJ "^^^ = [27) "^^^ = simpler ^^^^ 

In this notation, a = 1 means that the adiabatic and the isocurvature matter power spectra 
become the same in the long wavelength limit (see Eq.(|^)). 

As mentioned in the previous section, the isocurvature fluctuations give contributions to 
the COBE measurements which are larger than the adiabatic fluctuations by the factor /lo 
(if we take fig — 1 and h ~ 0.6, this factor is approximately ~ v^30)- Therefore, when we 
take account of the isocurvature fluctuations, the correct COBE normalization is expressed 
as ilTOll 



y3/2 



5.3 X 10- 



Ar=60 V 1 + /lo" 



Here, we have ignored the tensor perturbations. 
From Eqs. (O) and ([T8|), we see that 



^Our Eq.(l^) is different from Eq.(5) in Ref. 1) by a factor 1/4 due to a typo. 

^ The factor (2/15) comes from the value of the transfer function in the long wavelength Hmit, |2C] 



and the extra factor (10/9) is due to the decay of the gravitational potential at the transition from 
the radiation dominated universe to the matter dominated universe. 
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Thus, if we know the value of ot and if, we have a constraint on Fa- 



IV. INFLATION MODELS 

In order to determine the size of isocurvature fluctuations produced during inflation, here 
we consider some inflation models. 

There are mainly three classes of realistic inflation models: hybrid inflation models, new 
inflation models, and chaotic inflation models. We adopt simple realizations of the former 
two types of models in the context of supergravity. Although it is difficult to achieve chaotic 
inflation in supergravity, we also consider the chaotic type for comparison, since it is very 
simple and is often considered in the literature. 

In all cases, we have a reasonable parameter region which produces observable isocurva- 
ture fluctuations. 



A. Hybrid inflation model 

In this subsection, we consider the hybrid inflation model proposed by Linde and Ri- 
otto, 1^ because this hybrid inflation takes place under natural initial conditions and 
is consistent with supergravity (SUGRA) p2| (see also Refs. 23) and 24)). Let us now 
consider the hybrid inflation model, |^ which contains two kinds of superfields: S'(x, &) 



and '?/'(x, ff) together with '?/'(x, ff). Here Q is the Grassmann number denoting super- 
space. [|2^ The model is based on f/(l)_R symmetry under which S[ff) e^*"S'(6'e~*'^) and 
tlj{6)il!{6) ^/'(^e~*")'?/'(6'e~*"). The superpotential is then given by 

W = S{-fi'^ + Kipilj), (20) 

where /i is a mass scale and n a coupling constant. The scalar potential obtained from this 
superpotential, in the global SUSY limit, is 

12 



V 



+ D— terms, 



(21) 



where scalar components of the superfields are denoted by the same symbols as the corre- 
sponding superfields. The potential minimum, 

,2 



(5) = o, ww = ^, \m = m 



(22) 



lies in the D-fiat direction {ipl = By the appropriate gauge and /^-transformations in 
this D-fiat direction, we can bring the complex S, ip and ip fields onto the real axis: 



S = —;=cr, ih = ih = —c 



(23) 



where a and are canonically normalized real scalar fields. The potential in the D-flat 
directions then becomes 



1 \2 1 

2 I j,2 \ I -"^ 2 2i2 



(24) 



and the absolute potential minimum appears at a = 0, = = yu/y^- However, for 
0" > (Jc = \/2^/^/k, the potential has a minimum at = 0. The potential Eq.(^) for 
= is exactly flat in the a-direction. The one-loop corrected effective potential (along the 
inflationary trajectory a > ac with = 0) is given by p 



one— loop 



1287r2 



[K,a — 2/i ) m ■ 



A2 



2, Ka'^ + 2fi'^ 



+ {na^ + 2/^2)2 In 7. ^ - 2«:V^ In 



(25) 



A2 A2 

where A indicates the renormalization scale. 

Next, let us consider the supergravity (SUGRA) effects on the scalar potential (ignoring 
the one-loop corrections calculated above). The i?- invariant Kahler potential is given by ||26 



4 



(26) 



^We have assumed a U{1) gauge symmetry, where il){x^ 6) and iIj{x, 9) have opposite charges of 
the U{\), so that the il^ip term is allowed. 
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where the eUipsis denotes higher order terms, which we neglect in the present analysis. Then, 
the scalar potential becomes [E3] 



K a 1 



4 K 



2 4^ 

As in the global SUSY case, for a ^ (ic the potential has a minimum at 
potential for a ^ ac and = becomes 

VsuGRA = 1 ^ + 2 16 ^ J ■ 



(27) 

0. The scalar 
(28) 



In the first approximation, we assume that the infiaton potential for a ^ cXc and = 
is given by the simple sum of the one-loop corrections Eg. (pSj) and the SUGRA potential 
Eq.(|8D: 



V{a) = /iM 1 + ^(x' + 



/3 2 4/32 + 7/3 + 2 



16 



-a 



+ 



1287r2 



/ 2 o 2\2i '^'^^ 2/^2 

[Ka — 2fi ) in- 



A2 



+ («:a2 + 2/iyin^^^J-^-2«:Vln 

A2 A2 



(29) 



Hereafter, we study the dynamics of the hybrid inflation with this potential. 

We suppose that the infiaton a is chaotically distributed in space at the Planck time and 
it happens in some region in space that a is approximately equal to the gravitational scale 
and is very small (~ 0). Then, the infiaton a rolls slowly down the potential, and this 
region inflates and dominates the universe eventually. During the inflation, the potential 
is almost constant, and the Hubble parameter is given by if = 1^1/2/^3 ^ ^y^. When 
a reaches the critical value ac, the phase transition takes place and the inflation ends. In 



order to solve the flatness and horizon problem we need an e-folding number ^ 60. ||27 
In addition, the adiabatic density fluctuations during the inflation should account for the 
observation by COBE, which leads to Eq. (pTSp. 
The evolutions for a and A^ are described by 
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N = -H, (31) 

which are numerically integrated with cr|jv=o = c"c and Eq.(|TBp as boundary conditions. 
Though the inflaton potential is parameterized by three parameters (i.e., k, (3, and fi), we 
can reduce the number of free parameters from three to two by using the constraint Eg. (P^). 
We consider function of k and (3, i.e. /i = fi{K, (3). 

When cr ^ 1, the slow roll approximation cannot be maintained. Therefore, if the 
obtained value of o"|Ar=6o = o'o is larger than the gravitational scale, we should discard those 
parameter regions. Also, one of attractions of the hybrid inflation model is that one does 
not have to invoke extremely small coupling constants. Thus we assume all of the coupling 
constants n and (3 to have values of 0{1). Here, we choose 10^^ ^ i^, P ^ 10^^ as a 
"reasonable parameter region" (when k, /5 > 10^^, ctq exceeds unity, and cannot be as 
large as 60). 

The result is that H is as large as ~ (9(10^^^^^GeV), with the reasonable values of 
the coupling constants 10^^ ^ n, (3 ^ 10^^ (see Fig.^). When H has such a large value, 
the inflation generates large isocurvature fluctuations of the axion, if it exists. 

We are at the point to evaluate a, which depends on two parameters FaOa and H, as seen 
from Eg . ([19|) . We take H ~ lO^^^^^GeV, which has been obtained for the case of a = 0. 
We have, however, found that similar values of the Hubble constant H are obtained even 
for a 7^ as long as |a| ^ 1. From Eg.(0) we derive Fa9a ^ lO^^GeV for a ^ 0.1 and 
H ~ lO^^^^^GeV. Reversely, if a takes a value of (9(10^^) (for example, a ~ 0.05) as we 
have seen in §|I|, and H ~ lO^^GeV, we have FaOa — 1-5 x lO^^GeV. 

In the standard cosmology. Fa has the upper limit Fa ^ lO^^GeV. However, as shown 
in Ref. (28) this constraint is greatly relaxed if late-time entropy production takes place. In 
this case, the unclosure condition for the present universe leads to an upper bound for FaOa 
as 

FA ^ 4.4 X lO^^GeV, (32) 
15 



where the reheating temperature after late-time entropy production is taken as Tr = IMeV. 
This relaxed constraint is consistent with Fa (— 10^^ GeV) required for a ~ 0.05 and 



B. New inflation model 

In this subsection, we consider the new inflation model proposed by Izawa and 
Yanagida, which is based on an R symmetry in supergravity. 

In this model, the inflaton superfield (j){x,6) is assumed to have an R charge 2/{n + 1) 
so that the following tree-level superpotential is allowed: 

W = (33) 

n+1 ^ ^ 

where n is a positive integer and g denotes a coupling constant of order 1. We further assume 
that the continuous U{1)r symmetry is dynamically broken down to a discrete Z2nR at a 
scale V, generating an effective superpotential: 

W^jf = v'^-^r^\ (34) 

The /^-invariant effective Kahler potential is given by 

^(0) = i0r+^i0r+---, (35) 

where C is a constant of order 1. As shown in Ref. 29), the spectral index of the density 
fluctuations is given by 



1 - 2C. (36) 



By using the experimental constraint ^ 0.8, we take a relatively small value for C? C ~ O-l- 
Let us now discuss the inflationary dynamics of the new inflation model. We identify the 
inflaton field (f{x)/^/2{> 0) with the real part of the scalar component field of superfield (j). 
(We use the same symbol for a scalar component field as the superfield.) The potential for 
the inflaton is given by 
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,4 >>„A,„2 9 „,2,„n I 9 ,„2n 



^ t;^ - I^V^ - ^^^V" + (37) 

for < (^) = V2(0) = v^(i;V^)^/". Here, ^ and V are taken to be positive. It is shown in 
Ref. 29) that the slow-roll condition for the inflaton is satisfied for < 1 and if ^ (pf, where 



1 



This provides the value of at the end of infiation. Hereafter we take n = 4, since it is 
shown in Ref. 29) that this is the most plausible case. 

We assume that the infiaton begins rolling down its potential from near the origin ^ 
0). The Hubble parameter during the inflation {0 < ip ^ (ff) remains almost constant and is 
given hy H = V^^'^/y/S ^ v'^/^/S. When Lp reaches (pf, the inflation ends. As in the hybrid 
case, we need an e-folding number ~ 60, and we solve the equation of motion numerically. 
Also, we can flx one parameter v by using Eg. (p!8|) . The result is that H is as large as 
H ~ (9(10^^'')GeV with the reasonable parameter region 10^^ ^ g ^ 1, 10^'^ ^ C ~ 10~^ 
(see Fig. pup. This value is for a = 0, but similar values are obtained for a ^ 0.1. If we take 
a ~ C(10-2) and H ~ O(10^-^)GeV, we have from Eq.(^) F^Oa ~ lO^^GeV, which is the 
canonical value for axionic dark matter without late-time entropy production.^ 

C. Chaotic inflation model 

It is difficult to construct a realistic model of chaotic inflation in supergravity. However, 
a monomial type of chaotic infiation model is very simple and is treated widely in the 
literature. We consider such a model in this subsection. 

For simple models of chaotic infiation such as 



^For large C 0.05), the spectral index deviates from 1 (see Eq.(36)), and the observational 
constraint on a discussed in §|| becomes slightly stringent. However, the result in §0 can be 
directly applied for smaller (. 
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V{ct>) = ^ (39) 



or 



^(0) = (40) 

the COBE normalized value of Hubble constant during inflation is about H ~ (9(10^^)GeV. 
In this case, a too great amount of isocurvature fluctuations is produced, even considering 
late-time entropy production. However, if the potential for Peccei-Quinn scalar fleld is ex- 
tremely flat, the effective value of Fa during inflation can be larger and Fa ~ lO^^GeV is 
possible. [0 When this is the case, we have an interesting amount of isocurvature fluctua- 



tions again. For example, if Fa ~ 2.4 x lO^^GeV and a ~ 0.05, we have if ~ 1.4 x lO^^GeV, 
which is a natural value in chaotic inflation. fpJl 



V. CONCLUSIONS AND DISCUSSION 

In this paper we have studied density fluctuations which have both adiabatic and isocur- 
vature modes and have discussed their effects on the large scale structure of the universe and 
CMB anisotropies. By comparing the observations of the large scale structure, we have found 
that the mixed fluctuation model is consistent with observations if the ratio a of the power 
spectrum of isocurvature fluctuations to that of adiabatic fluctuations is less than ~ 0.1. 
In particular, the mixed fluctuation model with a ~ 0.05, total matter density Qq = 0.4, 
vacuum energy density Aq = 0.6, and Hubble parameter Hq = 70km/s/Mpc gives a very 
good flt to both the observations of the abundance of clusters and the shape parameter. 
(In recent observations of high-redshift supernovae, a similarly large value for Aq has been 
obtained. ||3^) Therefore, the mixture model of isocurvature and adiabatic fluctuations is 
astrophysically interesting. 

The CMB anisotropies induced by the isocurvature fluctuations can be distinguished 
from those produced by pure adiabatic fluctuations because the shapes of the angular power 
spectrum of CMB anisotropies are quite different from each other on small angular scales. 
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The most significant effect of the mixture of isocurvature fluctuations is that the acoustic 
peak in the angular power spectrum decreases. From observations of the CMB anisotropics 
in future satellite experiments, we may know to what degree the axionic isocurvature fluc- 
tuations are present. 

The requirement a ^ 0.1 leads to constraints on the Peccei-Quinn scale Fa and the Hub- 
ble parameter during the inflation. Isocurvature fluctuations with a ~ 0.05 are produced 
naturally in the hybrid inflation model if we take the Peccei-Quinn scale Fa — 10^^~^^ GeV. 
This Fa exceeds the usual upper bound {Fa ^ 10^^ GeV). However, this bound can be 
relaxed if there exists late-time entropy production, and in this case Fa ~ loi^"!^ GeV is 
allowed. Furthermore, Fa — 10^^ GeV is predicted in the M-theory. [|3^ Thus the exis- 
tence of large isocurvature fluctuations may provide crucial support for the M-theory axion 



hypothesis. 34 



In the new inflation model, we also have sufficiently large isocurvature fluctuations to 
be observed. In this case, we do not have to invoke late-time entropy production to dilute 
the axion density. In a simple chaotic inflation model, we have a larger value for H during 
inflation than the above two cases. In this case. Fa exceeds the upper limit, even considering 
the late-time entropy production. However, the Peccei-Quinn scalar field may have an 
extremely flat potential, and Fa may be as large as the gravitational scale. In this case, we 
have an appropriate amount of isocurvature fluctuations again. 

In all of these three realistic classes of inflation models, we have astrophysically interesting 
amounts of isocurvature fluctuations. 
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FIGURES 




IQl 102 IQ3 

I 

FIG. 1. A sample of the CMB anisotropies angular power spectrum (upper panel) and fi (lower 
panel). In the upper panel, the solid line corresponds to pure adiabatic fluctuations, and the short 
dashed line to pure isocurvature fluctuations. We have chosen Clo = 1, Ao = 0, h = 0.5, n = 1 
and nsh? = 0.015 here. 
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FIG. 2. as{ct) as functions of and h normalized by COBE. The four panels correspond to 
a = (top left), 0.02 (top right), 0.04 (bottom left) and 0.06 (bottom right). The hghtly shaded 
region inside the two short dash - long dash lines indicates the values of ag observed from the X-ray 
clusters of the galaxies. The three long dashed lines represent shape parameters T = 0.2, 0.25 and 
0.3 from left to right, and the observationally favored region (T = 0.25 it 0.05) is also lightly shaded. 
The darkly shaded region satisfies the constraints from both the cluster abundances and the shape 
parameter. 



23 



0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 
^^0 Q.0 



FIG. 3. The same as Fig. ^, except that the four panels correspond to a = 0.08 (top left), 0.10 
(top right), 0.12 (bottom left) and 0.14 (bottom right). 
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FIG. 4. The power spectra for a = 0.05, JIq = 0.4 and h = 0.7 (left panel) and a = 0,O,q = 0.3 
and h = 0.7 (right panel). The symbols denote the observational data, pl 
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FIG. 5. The solid lines represent values of a which give the best fit to the data of a% from 
cluster abundances. The three long dashed lines represent shape parameters V = 0.2, 0.25 and 0.3 
from left to right. 
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FIG. 6. Peak-height parameter 7 (see Eq.(|n])) for a = (top left), 0.04 (top right), 0.08 
(bottom left) and 0.12 (bottom right). 
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FIG. 7. The CMB angular power spectra normalized by COBE. Here we have chosen 
Qq = 0.4, h = 0.7, Ao = 0.6, and rj^/i^ = 0.015. The short dashed line corresponds to the 
pure adiabatic case (a = 0), the long dashed line to the pure isocurvature case {a = oo), and the 
solid line to the case a = 0.05, which gives the best fit to the observational data. 
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FIG. 8. The same as Fig. except that the cosmological parameters and h have been 
changed: VLq = 0.3, h = 0.9 (top left), Q.q = 0.4, h = 0.7 (top right), VLq = 0.5, h = 0.6 (bottom 
left) and VIq = 0.6, h = 0.5 (bottom right). In this figure, we also plot angular power spectra for 
various a to see how much the acoustic peak decrease as a changes. 
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FIG. 9. The Hubble parameter during the hybrid inflation, which is normahzed by the COBE, 
ignoring tensor perturbations and isocurvature fluctuations. In the region above the dashed hne, 
(To exceeds the gravitational scale (i.e., ctq > 1) and it is excluded. 
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FIG. 10. The Hubble parameter during the new inflation, which is normahzed by the COBE, 
ignoring tensor perturbations and isocurvature fluctuations. 
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